#ifndef UTILITI_H
#define UTILITI_H

//the thickness of the plate
#define H 0.001
//Young's modulus
#define E 2e11
//Poisson Ratio
#define GAMA 0.3

#include <iostream>
#include "Eigen/Core"
#include "Eigen/Eigen"

using namespace std;

Eigen::MatrixXd calculateD();

Eigen::MatrixXd calculateB(double x1, double y1, double x2, double y2, double x3, double y3, double x4, double y4, double s, double t);

double JacobiMatrix(double x1, double y1, double x2, double y2, double x3, double y3, double x4, double y4, double s, double t);

Eigen::MatrixXd calculateK(double x1, double y1, double x2, double y2, double x3, double y3, double x4, double y4);

#endif